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Abstract 



In this thesis we study dynamic of magnetic monopoles from Lagrangian density 
in Yang-Mills-Higgs field theory. In particular, we discuss BPS (Bogomolny Prasad 
Sommerfield) magnetic monopoles, described by SU(N) Bogomolny equations, which 
has field equations in form of non-linear coupled matrix field equations. One of 
the methods to simplify SU(N) Bogomolny equations is by using harmonic maps 
CP^ N ~ l \ This method has relation with Gr(n,N) cr-model and can transform SU(N) 
Bogomolny equation into more simple scalar field equations that depends only on one 
variable. As an example, we consider the case of SU(2) Bogomolny equation. 
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Chapter 1 
Introduction 



Magnetic monopoles concept was first introduced by P.A.M. Dirac while he tried 
to search for explanation about the unit of electronic charge e pQ. His research 
was based on the fact that electric charge is always observed in integral multiples 
of the electronic charge e. This electronic charge e has made % approximately 137 
which then became his focus of research and some others physicist, one of them was 
A.S. Eddington. In his research development, Dirac considered some arguments but 
unfortunately did not lead him to any value of 137 and for that reason he felt that 
his arguments were a failure. Instead, the result of his research, which he wrote it in 
a paper titled by Quantized Singularities in The Electromagnetic Field, born the new 
idea of magnetic monopoles. It was the concept that for years later interest many 
scientists with capable of wide generalizations. 

In the following years, many papers have been published in the topic of magnetic 
monopoles. Between year 1973 through 1976, there were more than 300 research 
papers on the subject of magnetic monopoles [2]. A primary contribution for the the- 
oretical investigations was provided by 't Hooft and Polyakov in 1974 that they dis- 
covered a magnetic monopoles solution in a spontaneously broken non-Abelian gauge 
theory [3] . Their work pointed out the natural manner in which magnetic monopoles 
make their appearance in these theories and encouraged further exploration of this 
phenomenon. 

Many studies about magnetic monopoles in non-Abelian gauge theory were moti- 
vated by two sources. First, if magnetic monopoles solitons occurs in a unified gauge 
theory of the weak and electromagnetic interactions, which then the theory describes 
the real world proved to be correct, then they will be experimentally accessible even 
though at extremely high energies. In this case, there must be available informa- 
tion about these soliton of magnetic monopoles field configurations. Second, there 
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may be a connection between magnetic monopoles and the quark confinement mech- 
anism. These magnetic monopoles would probably be counter parts of some charge 
other than electric charge and therefore be expected to have little connection with 
the magnetic monopoles sought in the experiments [3]. 

In this thesis, we describe how to develop SU(N) Bogomolny equations for BPS 
magnetic monopoles from Yang-Mills-Higgs field theory and applying a mathematical 
device, harmonic maps CP^ N ~ l \ to transform the equations in the more simple form 
that can be used for further methods to find the solutions. In the first chapter, we dis- 
cuss the Yang-Mills-Higgs field in Riemann sphere and derive its dynamics equations. 
Then by using its energy equation, we take some conditions and Bogomolny analysis 
to get SU(N) Bogomolny equations for BPS magnetic monopoles. For Chapter two, 
we investigate on harmonic maps and how its connection with Gr(n,N) cr-model. We 
take special case for CP^ N ~ l > space of Gr(n,N) cr-model and derive many of its prop- 
erties and awe also use a Veronese map as we will used for next chapter to obtain 
the solutions. Chapter three is the major work of this thesis where we use ansatz for 
SU(N) Bogomolny equations in chapter one to connect with harmonic maps method 
in chapter two. With that ansatz, we can transform the SU(N) Bogomolny equa- 
tions form its matrix form into scalar equations which are much more simple to be 
discussed. We also take an example for SU(2) Bogomolny equations and find its 
solutions. The last chapter is the conclusion about the results that we get in the 
previous chapters. 
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Chapter 2 

Magnetic Monopoles 



2.1 Dynamics of Yang- Mills- Higgs Field 

In searching for quantization of a unit of electric charge e, Dirac start his attempt by 
writing Maxwell equations in matter (with electric charge and current source) into 
symmetrical form between electric and magnetic field. This is based of the fact that 
in vacuum conditions, Maxwell equations has symmetrical form between its fields. 
So, why it does not occur in the same way as we move into Maxwell equations in 
matter. For that reason, Dirac introduces a magnetic charge and current source in 
his version of Maxwell equations pp. Next, this idea was extended in more general 
because as we know that Maxwell equations in vacuum basically can be derived from 
pure U(l) Yang-Mills-Higgs field theory. 

In this section, we discuss about magnetic monopoles which is derived from Yang- 
Mills-Higgs field theory with Lagrangian density given by [5] 

C = ~tr (F liV F' iv ) - Ur (D^D^) - ~ (||$|| 2 - l) 2 (2.1) 

where is the gauge field and $ is a Higgs field, as sources for magnetic monopoles, 
with F^y = d^Ay — d u A^ + and covariant derivative = + [A^,]. We 

also write the Higgs field in $ = i$ a T a which is a skew Hermitian form so that 
||$|| 2 = -\tr ($$) and with tr (T a T b ) = 25 ab (T a is generator of Lie group SU(N)). 
For this chapter, we use index convention //, v = 1, 2, 3 and a, b = 1, ... , dim[SU (N)]. 
If we look at the equation (2.1), we recognize than on the first part of equation (2.1) 
is for electric and magnetic field while the second part is for the sources of magnetic 
monopoles and then the third part is a Lagrange multiplier that comes from constraint 
condition ||$|| =1. In order to write dynamics equations from Lagrangian density 
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above, we have to take variation of action to each of its fields. As variation of action 
is wrote by 

5S = J 5C d A x (2.2) 

If we look for dynamics equations of field then we take variation of to field to 
the equation (2.2) and action on Lagrangian becomes 

SC = hr(d,(F^5A u )-d u (F^5A^-d,F^5A u + d u F^5A,-F^[A v ,5A^] 

+ F>" [Ap, 5A U }) + Ur [<t>, 5A,\) (2.3) 

Because we have boundary condition for integral on surface SA^ = 0, then the second 
part of equation (2.3) can be ignored so that 

5C = hr((D u F^ + [D»$, $])<L4 M ) (2.4) 

from least action principal SS = 0, then dynamics equation for fields A^ are 

D V F^ V + [£>"$, $] = (2.5) 

While for $ field, we take variation of $ field to the equation (2.2) and the action on 
Lagrangian becomes 

5C = --trid^D^5<$>) - d^<$>5<5> + (D^A^S® - D^SQA^) 

+ ^(||$|| 2 - l)tr($5$) (2.6) 

and also for integral on surface we have boundary condition <5$ = 0, then the second 
part of equation (2.6) can be ignored so that 

5C = hr((D,D^ + ^(||$|| 2 - (2.7) 

From least action principal 8S — 0, then dynamics equation of $ field is 

D^D^+ -(||$|| 2 - 1)$) = (2.8) 

As we can see in equation (2.5) which is similar to Maxwell equations in matter with 
current sources is in form of field $. 
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2.2 Energy Equation for Yang-Mills-Higgs Field 

Now, let we find energy equation for Lagrangian (2.1) and then use it to get SU(N) 
Bogomolny equations. To do that, we have to count the tensor energy- momentum for 
Yang-Mill-Higgs field by rewriting the action of Lagrangian (2.1) for arbitrary metric 

S = J d 4 x (2.9) 

where g = det(g^u) = e tr< - ln 9 ^"\ Then, we take variation over g pu to the action (2.9) 
so that 

5S = j{5{^=g~)C + ^5C)d A x (2.10) 

with 

SV = g='^{-9)~H-Sg) 
Sg = e^ ln ^hri—5g, u ) = g tr{g^5g» v ) (2.11) 
If we write in its components 5g = g g pu Sg pu , so that 

Sy/=9=l(-9)^6 gilv (2.12) 

While 

5C = ^tr(F pa F^)Sg^ - hr(D p $D^)5g^ (2.13) 
and using Sg^ p = —g m bg a p9^ p -, then 

5C = {- X -g»°tr{F pa F pv ) + hr(D^D^))g^g^5g^ (2.14) 
Substitute equations (2.12) and (2.14) into equation (2.10), then 
SS = -\^F~g J {-g aP C + l -{g™tr{F pa F, v ) - tr{D^D^))g^g^)5g aP d A x (2.15) 

From general relativity for action of matter is 5 S = —\yf—g J T a ^bg a ^x so we have 
tensor energy-momentum 

T a,3 = _ ga p c + 1 ( g ™ tr ( FpaFfiu ) _ tr(D p <S>D^))g» a gPP (2.16) 

and take its energy part T 00 = -g 00 C + ±(g" a tr(F pa F pu ) -tr(D p <$>D p <$>))g^g op . From 
definition, we know that the energy is written by E — J T 00 ^gd 3 x for arbitrary 
metric. Below, we write some metrics where ds 2 = g pu dx^dx u : 
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1. Standard Minkowskian metric 

ds 2 = dt 2 - dx 2 - dy 2 - dz 2 

2. Minkowskian metric with spherical coordinates 

ds 2 = dt 2 - dr 2 - r 2 d6 2 - r 2 sin 2 6d(f) 2 

3. Minkowskian metric with Riemann sphere coordinates 

2r 2 



(2.17) 
(2.18) 



ds = dt — dr 



9fj.u 



(1 + KI 2 ) 2 
/ 1 
0-1 


V o o 



(« + «) 






2r 2 



(l+l?l 2 ) 2 

derivation of metric (2.19) is in Appendix A. 
So, we have the energy part T 00 with metric (2.19) is 

rpOO 



(i+l?l 2 ) 2 




(2.19) 



-C + ^g™tr(F 0a F 0l/ ) - hr(D $D $) 



(2.20) 



2.2.1 Static field and BPS condition 



For this thesis, we will only look for solutions of static fields and in BPS limit. In that 
case, we have for static field (<9 $ = 0, A = 0) and with BPS condition on A = 0, 
then equation (2.20) becomes 

T 00 = -^Pg^tr(F pa F^) + ^g ij tr(F 0i F 0j ) + ^g l Hr(D^D^) (2.21) 

where we use index i and j — 1,2,3 — r, and for index fj,, v, p and a = 0, 1, 2, 3 = 
t, r, ££. We can calculate equation (2.21) separately for each field. For the first one, 
we calculate the part that contains $ field which is written as below 



-gVtrip&DjQ) 



2\2 



-tr D r $D r § + 



(i + iei 2 ) 



While for the others part that bear field which are 



- -g^g^trXFM 
we can separate it again as below 



-g ij tr(F 0i F j) 



(2.22) 



(2.23) 
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• fI = o^ u ^O j p = )( t ^ 
-IgVtriFojFoi) 

-l9 ij tr(F 0i F 0j ) 

^ -yik g n tr(yFklF ..) 

Joining again all of its parts of the equation (2.21), then we have 

T 00 = -\tr (D r QD r $ + {1+ jf )2 D^D^ + ^Hr(F 0z F 0j ) - ^g lk g> l tr(F kl F ZJ ) 

(2.24) 

2.2.2 Magnetic monopoles conditions 

For magnetic monopoles it means that the sources only produces magnetic field or we 
may write F 0i = (if F 0i ^ then it is known as dyons,it means the sources produces 
two fields which are magnetic and electric field). With this condition then equation 
(2.24) becomes 

(2.25) 

So, we have equation of energy for BPS magnetic monopoles as below 

1 r ( 2r 2 (14- \f\ 2 ) 2 \ 

E = —Jtr ^ i + m2 D r *D r * + + 2F r ,F ri + L_J|LL % ^ j ^ 

(2.26) 

2.3 Magnetic Monopoles Equations of Motion 

From the energy (2.26), we can derive equations of motion for each field by using 
variation of energy about each field and least action principal SE = 0. So, for $ field, 
the variation of the energy (2.26) is 



2r 2 \ / 2r 2 



SE = ~Jtr(d r ( . .!> j - <),■ ( "' , , l> <\- \ H> 
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2r 



(2.27) 



(1 + lei 2 ) 2 

+<%( J D C <M$) - %(L> e $)<5$ + [D^,A^)d 3 x 
Because we have boundary condition for surface integral 5(f) = 0, then 

2r 2 



and equation of motion of $ field is 
/ 2r 2 



D, 



D r $J - Dz(Df$) - D^(D^) = 



v (l -K 

D r (r 2 D r <S>) = - (1+ ^' 2)2 (^(^) + Dz(D&j) 
With the same manner for A% field, then variation of the energy becomes 

5E = -1 /" tr( - [Df$, $\6A + d r (F r ^5A) - d r F ri 5A + [F rf , A r ]5A 



(2.28) 

(2.29) 
(2.30) 



(i+Jj! 

2r 2 



2\2 



(1 + KI 2 ) 2 
2r 2 



(2.31) 



and using boundary condition for surface integral 5 A — 0, then 

SE = -^J tr (-[D?f>, $}5A - D r F r ^5A + D ^ 1 ^ F^)5Aj d 3 x (2.32) 
So, equation of motion for A^ field is 

- [D&, $] - D r F rf + £> f 



(1 + ld 2 ) 
2r 2 



^<t>, $] + D r F r? - = -L^- ((1 + |£| 2 ) 2 F ft ) 



(2.33) 
(2.34) 



If we do for field then we have the same equation of motion above just by changing 
the index £ -h- £ from equation (2.34), so the equation of motion for field A^ is 



[£> € <&, $ ] + A~F rf = ^ ((1 + |£| 2 ) 2 ^-) 



(2.35) 



While for A r field, if we take its variation on the energy, then 

,2 

+d^ r 5A) - d^ r 5A + [F fr , A e ]<L4) rf 3 x (2.36) 



5£ = -\jtr(- [D r $, $]5A + d^ r SA) - d^ r SA + [F €r , A { -]<L4 



and with boundary condition for surface integral 5A r , then equation (2.36) becomes 
8E = ~J tr (~ (1 ^| 2)2 [A-$, - £>^ r 5A - I^ r <S^ rf 3 a; (2.37) 



So, equation of motion for A r field is 

2r 2 



I + | , |2 )2 -[^' $ ] " ^ - ^fr = ( 2 - 38 ) 



[D r $, $] = (1+ Jf 2)2 (l> € -i^ + £> € i^-) (2.39) 



2.3.1 SU(N) Bogomolny equation 

In the previous, we have derive equations of motion (2.30), (2.34), (2.35), and (2.39) 
which are known as SU (N) non-Bogomolny equations for BPS magnetic monopoles. 
Those equations are in second order differential and have properties of non-linear and 
coupled fields, and it is hard to find solutions for that kind of equations. Next, we 
will derive SU(N) Bogomolny equations for BPS magnetic monopoles which have 
equations in form of first order differential. Those kind of equations are more simple 
then SU(N) non-Bogomolny equations and it is possible to find the solutions. In 
deriving SU (N) Bogomolny equations, we use Bogomolny analysis for BPS magnetic 
monopoles to the energy equation (2.26). In Bogomolny analysis, then energy equa- 
tion (2.26) is transformed into equation that contains square sum of the fields. We 
must consider that the energy has to have real value E > 0, so that : 



1. For first part of the energy 



= -<9 r $ + [4, $] (2.40) 



then At = -A r . 
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2. and second part 



(d& - [4 *]) (2.41) 



then = -A^. 



3. and third part is correct by using the result from 1) and 2) 



- F r £ — [F^ 



= -{d r A^-d^A r + [A r ,A^}) (2.42) 
4. also for fourth part is correct if we use the result from 1) and 2) 



= + (2.43) 

Let we look for some form of equations below: 

(iD^ - F r t)(iDz$ - F ri y = -D^D^ - F ri F ri + iD^F^ - iF^D^ (2.44) 

then we can transform energy equation (2.26) so that it contains equations (2.44) and 
(2.45), and we may write it as 

|2\2 



i f ,(i + \£\ 2 y 

E = lj tr V r 2 (iD& - F r e)(iD& + F ri ) 



( , 46) 
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This energy (2.46) has minimum a value if 

iD^ = F ri (2.47) 
iD r S> = 1 } % (2.48) 

and those equations are known as SU(N) Bogomolny equations for BPS magnetic 
monopoles in Riemann sphere metric (2.19). 
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Chapter 3 



Harmonic Maps 



3.1 Harmonic Maps Definition 



Basically, harmonic map is a map between Riemannian manifolds. The theory of 
harmonic maps was introduced in 1945 by B.F. Fuller and developed by J. Eeels and 
J. M. Sampson ten years later [7]. Its role in physics was shown firstly by C. W. 
Misner [8] in 1978, who formulated non-linear a-model field theory in geometrical 
description. The non-linear cr-models in two dimensions are special interest because 
they bear many similarities to the non-abelian gauge theory in four dimensions and 
have a property of being an integrable system. Later in this thesis, the non-linear 
cr-model is shortly named cr-model. 

Next, we explain some basic definitions of harmonic maps theory. Let Aio is a 
Riemannian manifold (source manifold) with local coordinates u % and metric 



with i, j = 1, 2, ... , dim[Aio] and M. is another Riemannian manifold (target mani- 
fold) with local coordinates f A with metric 



ds 2 = gijdtfdu-i 



(3.1) 



ds 2 = h AB df A df 



B 



(3.2) 



with A,B = 1,2,..., dim[M] 
Then, a map 



f:M ^M 

f = m 



(3.3) 
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is called a harmonic maps if it extremizes the action 



S= / d n °u^g~C (3.4) 

J Mo 



such that SfS = with n = dim[Mo] and g = det(gij) where 



1 df A df B •■ 

£ = ^ h,AB J ■ I - g v (3.5) 

2 <9w* ^ y v ; 



is the Lagrangian density. 



3.2 cr-Model 

The theory of harmonic maps in physics literature is related to cr-model. The a-model 
is a field theory with the following properties [9]: 

1. Fields 4>(u) of a-model are constraints for all points <p(u), Vm G A4q. 

2. The constrains and Lagrangian density are invariant under action of global 
symmetry group G on <j)(u). 

3.2.1 0(N) a-model 

As an example, we look on O(N) cr-model which consist of iV-real scalar fields (p A 
with A = 1, . . . , N and has Lagrangian density 

2 du l du^ y ' 

where scalar fields (p A satisfy the constraints 

<j> A (/) A = 1 (3.7) 

with repeated index means the sum of all its value. Lagrangian density (3.6) is 
invariant under global transformation group O(N) 

<p A -> (<j> A y = oU* (s.s) 
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Geometrically, the constraints (3.7) define a sphere (N — 1) dimension S N 1 in N 
dimension Euclidean space R N of fields manifold <p A . This constraints can be solved 
by introducing parametrization below: 



cf> A = f A , 7V = ± V / 1-|/| 2 , A = 1,2,...,(N-1) (3.9) 

where \ f\ 2 = f A f A and the range of |/| is limited to —1 < |/| < 1. The choice of sign 
+ or in equations (3.9) determines parametrization whether we use upper or lower 
hemisphere of S N-1 . By parametrization in (3.9), we can write Lagrangian density 
(3.6) as 

£=5^34* (3-10) 

with Hab is tensor metric of target manifold S^N — 1). The tensor metric Hab is 
defined by substituting parametrization (3.9) into metric 

d 2 a = d<fi A d(j) A 

= (*AB + -/4^) d f Ad f B (3-H) 



l/l 



so that 



h A B = 5 AB + -^- (3.12) 
1 - \f\ 2 

where 5ab is the Kronecker delta. As the result, the fields / are free of constraints 
and are solution of O(N) cr-model field, so its defines a harmonic maps 

f:M ^S N -\ (3.13) 



3.2.2 Group formulation in cr-model 

In group formulation, target manifold Ai of a cr-model, with G is the invariant global 
group, are homogeneous space of G. It means that by applying transformation group 
G over a field (f> p G M. then we can get into all the fields in target manifold M.. In 
another words, for an arbitrary field <p g ^ 4> p G M. at least there is an element g G G 
such that (j) q = gcf> p . 

If there is any stability or isotropy group H C G on the fields 4> p G M, 

H = {he G\h(f) p = <f) p } (3.14) 
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then target Manifold Ai is a coset space of G/H = {gh\g e G} that works on the 
field (f) p . If the identity Iq is the only subgroup of G, then Ai = G/H = G is manifold 
of group G, then we call it as chiral model. 

In general, the target manifold M, of cr-model is a manifold of group G. In that 
case, the cr-model can be represented in group G by writing parameters of group G 
as the fields f A , A = 1, . . . , (n — dim[G]), with dim[G] is dimension of group G or 
number of generator in group G. Hence, the tensor metric for target manifold in 
cr-model is written by 

**» = ~ 2tr ( G ~w G ~w) (315) 

and the Lagrangian density is 

£ = -* { G -w G 'w) (316) 

3.3 Grassmanian cr-Model 

Grassmannian manifold Gr(n,N), 1 < n < N, is the manifold of n dimensional 
planes passing through the origin in the iV dimensional complex space C N , where 
Gr(l, N) = CP N ~ l . So, the CP 1 ^" 1 is a set of lines passing through the origin in iV 
dimensional complex space C N such that a point in CP^" 1 is a line on C N . In this 
thesis, we use Grassmannian manifold in complex coordinates. 
In coset space formulation, the Grassmannian manifold is 

G ^= m^W) (3 - 17) 

and if we take into account the orientation of the planes, then 

^•^- SUiN-VsUin) < 3 - 18 > 

Gr(n,N) cr-Model consists of (N x n) complex matrix fields Z = (Z Aa ) with A = 
1, . . . , N and a = 1, . . . , n. Those fields satisfy the constraint: 

Z^Z = I n (3.19) 
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The Lagrangian density of the model is required to be invariant under global unitary 
transformations G G U (N) that work on Z field from the left 

Z -> Z' = GZ (3.20) 

and also invariant under local gauge transformations H(x) G U(n) from the right of 
Z field 

Z ->■ Z' = ZH{x) (3.21) 
In this thesis, we use Lagrangian density given by W. J. Zakrzewski [10] 



£ = tr ([D^Z^D^Z) (3.22) 

with 

D^Z = d^Z - ZZ ] d^Z (3.23) 

Next, we define w a , a = 1, . . . , n to be a set of N components orthonormal vectors 
in C , then the matrix field Z can be represented as 

Z = (w\...,w n ) (3.24) 

In this representation, Z defines an orthonormal n-frame in C . 
Let we choose Z as the last n-columns of matrix G 

G = (Y AB Z Aa ) , B = 1, . . . , (N - n) (3.25) 

where Y is an (N x (N — nj) matrix and Z is an (N x n) matrix. Since we have 
unitary condition for G G U(N) which is G^G = GG^ = In, then we have 

yty = Ilf _ nj y\ Z = 0, Z f Z = I n , YY^ + ZZ f = I N (3.26) 
that makes Z satisfies the constraint (3.19). For our purpose, we define field 

= {(Y Z) - (0 2Z)} ( Y ^\ (3.27) 

and by using the property in equation (3.24), then we can write as 

$ = (I - 2P) (3.28) 

where P = ZZ^ . A complete discussion of Gr(n, N) a- model and for interested reader 
can study the reference on dissertation of Hans Jacobus Wospakrik pTj . 
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3.4 Gr(n, N) cr-Model in Projection Space 

The Lagrangian density for the Gr(n,N) cr-model is simply given by W. J. Za- 
krzewski [10] 

C = -tr(d^d^) (3.29) 
8 

with $ = / — 2P and P is the projection operator, then 



£ = htr(d^Pd» P) (3.30) 



Taking into account the constraint below: 

P 2 = p = pt 

(P 2 -P) = (3.31) 
then we can write Lagrangian density as 

£ = hbr^Pdi* P) + \tr{P 2 - P) (3.32) 

So, action on Lagrangian density (3.32) is 

S = J (htr(d^Pd» P) + Xtr{P 2 - P)\ d^x (3.33) 

3.4.1 Equations of Motion for Gr(n,N) cr-model 

Equations of motion of Lagrangian density (3.32) is searched by variating the action 
(3.33) about P and A, then we get 

SS = Jtr (d^d^PSP) - d^d^P5P + A(2P - I)8P + (P 2 - P)5\) d 4 x (3.34) 

Because of boundary condition 5P = on surface integral, we can discard the first 
part of equation (3.34). And then, we get the equations of motion as below 

d^P - A(2P - /) = (3.35) 
P 2 - P = (3.36) 
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We can see that the equation (3.36) returns to the constraint (3.31). Next, if we 
multiply the equation (3.35) by P on the right and left separately, then we have two 
equations 

d^d^PP - A(2P -I)P = (3.37) 
Pd^d^P - XP(2P - I) = (3.38) 

If we count the difference between equation (3.37) and equation (3.38) then use prop- 
erty (3.31), then we get 

[P,& i d li P]=0 (3.39) 



3.4.2 Gr(n,N) cr-model solution 

In this section, we discuss the construction of solutions of Gr(n,N) a- model from 
equation (3.39) in two dimensional Euclidean space R 2 or the complex plane C. 
Later, we compactify R 2 by including points at oo to obtain the Riemann sphere 
S 2 = R 2 U oo and consider the harmonic maps: S 2 — > Gr(n, N). For next discussion, 
we use a Grin, N) cr-model that has source manifold in two dimensional Euclidean 
space. 

In searching for the solutions, we use the complex coordinates (£, £) for Gr(n, N) cr- 
model therefore the equation (3.39) becomes 



P, <%<%P 



= 



(3.40) 



Derivation of equation (3.40) is in Appendix B. We also can write equation (3.40) in 
form of 

[P, %P] + <% [P, <%P] = (3.41) 
Complete derivation of this equation in in Appendix B. 



3.4.3 Instanton solution 

Next, we look at differentiation of operator property P 2 = P on d^(P 2 = P) and 
<9e(P 2 = P) then we get 

d^PP = d ( P - Pd ( P (3.42) 

P<%P = %P - %PP (3.43) 

while the equation (3.41) is written by 

(P%P - <%PP) + % (P<%P - <%PP) = (3.44) 
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Substitute equations(3.42) and (3.43) to equation (3.44), then we get 

<% (Pd^P) - d 6 (d^PP) = (3.45) 

Special class of solutions for equation (3.45) is satisfied when Pd^P = such that 
(Pd^P)^ = d^PP = 0, so Pd^P = is the solution of field equation which is called 
selfdual equation. If we do the same way by substituting Pd^P and d^PP in equations 
(3.42) and (3.423) into equation (3.44), then it will produce a solution PdgP = which 
is called anti- selfdual equation. 

In this thesis, we focus on solution from selfdual equation Pd^P = 0. For this 
purpose, we consider an un-normalized (N x n) matrix field M = M(£, £) for which 
the \M\ 2 = M^M is assumed to be non-singular. As \M\ 2 is Hermitian, its eigen 
values are real and there exist an unitary matrix U such that 

|M| 2 = U*h 2 U (3.46) 

where A 2 is a diagonal matrix with eigen values (A 2 , . . . , X%). We define A to be the 
square root matrix of A 2 , hence A = (A 2 ) 5 = (Ai, . . . ,X N ). So, matrix |M| 2 = WA 2 U 
can be denoted by |M| 2 = (WAU) (WAU) such that \M\ = (|M| 2 )^ = WAU. 
In terms of M, the matrix field Z of Gr{n, N) cr-model is given by 

Z = M\M\~ 

and so the P projection operator is 

P = ZZ ] = (M\M\- 1 ) (JMI^MT) = M|M| 2 M f 

where we have use property of matrix (|M| -1 )^ = |M| _1 
Using property of 

(/ - P)M = (3.49) 

then the selfdual equation is given by 

M\M\- 2 (d^M)\l - P) = $ (3.50) 

Complete derivation of equations (3.49) and (3.50) is available in Appendix C. 

The solution of equation (3.50) is d^M = or M = M (£) which is called a holo- 
morphic matrix field. This solution of selfdual equation is named instanton solution. 
With the same analogy, if we discuss about anti- selfdual equation then we will get 
anti-holomorphic matrix field M = M (£) which is named anti-instanton solution. 
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(3.47) 



(3.48) 



3.5 Full Solutions of Gr(n, N) a- Model 



In this section, we will discuss the method of generating more general exact solutions 
of the two dimensional Gr(n,N) a- model starting from an instanton solution. The 
method was originally introduced by A. Din and W. J Zakrzewski [T2| [TO] . 

Let M k = M k (£, 0, k = 1, . . . , A where A < (N - l),be a set of (A + 1) mutually 
orthogonal (N x n) matrices, with n < N 

MjM, = \M k \ 2 5 kl (3.51) 

where 

\M k \ 2 = M\M k (3.52) 

are (n x n) non-singular matrices. Then the corresponding projection operator for 
each matrix M k is given by 

P k {n) =M k \M k \- 2 M\ (3.53) 

Clearly, tr(P k (n)) = tr(I n ), which means that each projection operator has rank- 
n. From equation (3.53), we can show that the projection operators are mutually 
orthogonal and also Hermitian 

P fc (n)P z (n) = 8 H Pi{n) (3.54) 
P k (ny = P k (n) (3.55) 

where \M\ 2 are Hermitian. In the following we want to present a generalized har- 
monic maps ansatz.To do this we use a sequence of mutually orthogonal matri- 
ces (M ,Mi, . . . ,M\) obtained from a sequence of holomorphic (analytic) matrices 
(M, <9 C M, . . . , <9|M), <%M = through the process of Gram-Schmidt orthogonaliza- 
tion. 

We can do this using the operator P + which is defined by its action on any matrix 

M e C Nxn H2UTDI 

M ikM MM 2 M ih.M (3.56) 
Then we have M M. .U, P+M, . . . , M k = P*M, . . . , M\ — P$M or simply 

M = M, M k = (/ - P^dzMk-u k = 1, . . . , A (3.57) 

where P k -i is the projection operator (3.53). An equivalent formulation for the se- 
quence M k , in terms of the projection operators P k , is given by 

M k = (1 - P Pk-i)d k e M (3.58) 
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With either one of these constructions in equations (3.57) and (3.58) the following 
properties of the matrices M k hold when Mq is holomorphic 

d- i M k = -M k _ 1 \M k _ 1 \- 2 \M k \ 2 (3.59) 
8 (M k \M k \~ 2 ) = M k+1 \M k \- 2 (3.60) 

Detail derivation is in Appendix D. 

Notice also that, for the CP^ N ^ l "> case, projection operators P k with k = 0, . . . , (N— 
1) are complete 

Po + Pi + ■ • ■ + Pn-i = I (3-61) 
and according to the construction M k in equation (3.58), then 

M N = (3.62) 

With the projection operators P k that we have constructed above, we have the 
following result that was originally proved by A. Din and W.J. Zakrzewski [12] using 
the Z fields formalism. 

Theorem: Each (NxN) projection operator P k {n) = M k \M k \~ 2 Ml, k = 0, 1, . . . , A, 
where M k = P£M and M = M (£) is a holomorphic (N x n) matrix field, solves 
the Gr(n,N) a-model (2.19) for its equation of motion (3.40). 

As a proof, we use the method that was developed by Sasaki [13J, where we shall 
be using the projection operator in term of M k matrix field as follows: 

From properties (3.58) and (3.59), we derive: 

d^Po = d ( {M Q \M \- 2 )Mi + M Q \M \- 2 d e Mt 

= MilMoT 2 ]^ (3.63) 
dzP k = d i {M k \M k \- 2 )Ml + M k \M k \- 2 {d-^M k y 

= M k+1 \M k \- 2 Ml- MfclAffc-il^MjL (3.64) 

with k = 1, . . . , X. 
Define 

fc-i 

Q k = J2 P h 2o = (3.65) 

1=0 

and take its differential about £ then use equations (3.63) and (3.64) 

fc-i 

1=0 

= M k \M k ^\- 2 Ml_ x (3.66) 
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Using orthogonality property of M k we obtain 

Q fe <9 5 Q fe = (3.67) 

where P k satisfy 

Q k dsP k = 0, P k d^Q k = dzQ k (3.68) 

Then define operator 

Rk = Qk + Pk (3.69) 

which satisfy self-dual equation 

R k d x iK = dsQk + Pkd^Pk 

= (3.70) 

and take the Hermitian conjugate of (3.70) gives 

d i Q k + d i P k P k = (3.71) 

If we take differential of equation (3.70) about £ and of equation (3.71) about £, 
then we get 

d&Qn + d^P k d^P k + P k d&P k = (3.72) 

and 

d^Q k + dzP k dtP k + d&P k P k = (3.73) 
The difference between equation (3.72) and (3.73) gives us 



P k ,d&P k =0 (3.74) 



which completes the proof. 



3.6 Action of Full Solutions and Non- Abelian Toda 
Equations 

As each projection operator P k (n) solves the Gr(n, N) a- model as stated in previous 
theorem (3.74), each projection operator P k (n) describes a specific field configuration 
having the action 

S k = % J « tr (%P fc <9 ? P fe ) (3.75) 
where So corresponds to instanton(or anti-instanton) configuration. 
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Using relation relation in equations (3.59) and (3.60), we have 
«%P fc = d-^M k \M k \~ 2 Ml + M k d-A\M k \~ 2 Ml 



= -M k ^\M k ^\- 2 Ml + M k \M k \- 2 Mt +1 (3.76) 

If we substitute equations (3.64) and (3.76) into equation (3.75), then the action 
becomes 

S k = % J <%d£ tr (|M fc | 2 |M fc _i|- 2 + |M fe+1 | 2 |M fe |- 2 ) 

= 2vr(A4+A4-i) (3.77) 

where 

Mh = hS d ^ tr (\ M ^\ 2 \ M ^ 2 ) ( 3 - 78 ) 

and by definition ]\f-\ = because d^M = 0. 

In the following, we derive recurrence relations for tr (\M k \ 2 \M k _i\~ 2 ) appearing 
in the integral (3.78). To do this we rewrite the definition of M k+1 in (3.57) as below: 

dzM k = M k+1 + P k d^M k (3.79) 

By using d^M k in equation (3.59) and property 

M^M k = -MlM^M^-^Mkl 2 

= (3.80) 

then we take differential of (3.79) about £ and multiply on the left by M k and on the 
right by |M fe |~ 2 , and becomes 

Mld&M k \M k \~ 2 = Mld-^M k+l \M k \~ 2 + Mld-^P k d^M k ) \M k \~ 2 (3.81) 

The left side of equation (3.81) gives 

Mld&M k \M k \- 2 = -\M k \ 2 \M k _,\- 2 (3.82) 

.From equation (3.59), the first part of the right side gives 

Mld^M k+1 \M k \~ 2 = -M\M k \M k \- 2 \M k+l \ 2 \M k \- 2 

= -\M k+l \ 2 \M k \- 2 (3.83) 
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and the second part becomes 

Mld^{P k d^M k ) \M k \~ 2 = d i (%|M fc | 2 |M fc |- 2 ) (3.84) 

The detail of equations (3.82) and (3.84) is in Appendix E. 

Substitute equations (3.82), (3,83), and (3.84) into equation (3.81) that gives 

(d^\M k \ 2 \M k \- 2 ) = \M k+1 \ 2 \M k \- 2 - \M k \ 2 \M k _i\~ 2 (3.85) 

For the n = 1 or CP(N-l) case, when M k are s sequence of iV-component vector fields, 
the equation (3.85) gives the Toda equation [H]. If we take the trace of (3.85) then 
we obtain 

tr (d ( (d^\M k \ 2 \M k \- 2 )) = d&tr (ln\M k \ 2 ) (3.86) 
Because det\M k \ 2 = e tr ( ln \ Mk \ 0, we may write equation (3.86) as 

tr (<%|M fc | 2 |M fc |~ 2 )) = d^ln (det\M k \ 2 ) (3.87) 
So, the trace of equation (3.85) is 

d&ln (det\M k \ 2 ) = tr (|M fc+1 | 2 |M fe |- 2 ) - tr (jM fc | 2 |M fc _i|~ 2 ) (3.88) 
Note that, for Gr(l,N) = CP N ~ l case, the equation (3.79) becomes 

d&ln (\M k \ 2 ) = \M k+l \ 2 \M k \~ 2 - \M k \ 2 \M k ^\' 2 (3.89) 

3.7 Veronese Map 

In this subsection, we describe how to construct the full solutions of the two di- 
mensional Grassmanian a-model using Veronese map. Here, we consider only for 
Gr(l,N) = CP N ~ X case with Mq G C Nx1 is a N-component vector field in complex 
space 

M : C ->• C N 

Z^M = (f ,...,f p ,...,f N ^) T (3.90) 
In the Veronese map [15] . we choose f p as function of £ in order p such that 

|M | 2 = (l + lei 2 )"" 1 (3.91) 
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This condition restricts the component of f p to have the form f p = ^JC^ where 
C^ -1 is the combinatorial factor and so the corresponding CP^ -1 field is 



M _ (' g 



, T 
iV-l x 



iMn 



1 + 



(3.92) 



For general case N in the sequence of M&, we use proof that was first given by 
Ioannidou et. al [16] : 

Proposition:For the Veronese map Mq in (3.92), the sequence Mk, (k = 0, . . . , N— 
1), constructed by equations (3.57) or (3.58), called Veronese sequence [15J, satisfy 

|M fc+1 | 2 (k + l)(N-k-1) 



\Mk\ 2 (1 + KI 2 ) 

To prove it, we make use of the recurrence relations (3.85) or (3.88). As |M | 2 = 
(1 + \m N ~\ then 

d^ln\M \ 2 = (N - 1) (l + |£| 2 )~ 2 (3.94) 
and by using relation (3.89) for k = 0, we obtain 



d^ln\M \' z = |Mi| 2 |M r^ - |M HM_! 



I Mi 



|2 



|M | 2 

where by definition M_i = 0. From equations (3.94) and (3.95), we have 

(Mil 2 _ (JV-1) 

PoF~ (i + \m 2 

For general fc, 1 < k < (N — 1), we use inductive proof by assuming 

|M fc | 2 k(N-k) 
l^-il 2 ~ (1 + KI 2 ) 2 
which is already true for k = 1. We can write |M&| 2 using equation (3.97) as 



(3.95) 



(3.96) 



(3.97) 



„ IMJ 2 \M k A z Ml, . , r 



2 

2 



M fc _ 1 | 2 |M fc _ 2 p'"|Mo| 21 
fc!(iV-l)! / „„\M-i 



(N-k- 
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then if we take differential over £ and £ of equation (3.98), we obtain 

d^ln\M \ 



(N-2k- 1) 
' (1 + KI 2 ) 2 



(3.99) 



Substitute equations (3.97) and (3.99) into (3.89), gives 

(N — 2k — 1) |M fc+ i| 2 k(N-k) 



a + iei 2 r 



\M k 



a + iei 2 r 



therefore 



\M k+1 \ 2 (k + l)(N-k-l) 



(3.100) 



(3.101) 



(1 + Kl 2 ) 2 
is true. 

Corollary: Each configuration Mk, (k = 0, ...,N — 1), of the Veronese map 
(3.92), has 

Af k = (k + l)(N-k-l) (3.102) 
We prove it using the integral formula on surface 



fa 



a + iei 2 r 



% \ det 



e i(t> ipe** 
e -i * -ipe i4 



dpdcj) 
[l+p 2 ) 2 



= 2ir 



(3.103) 



with p = tan (I). Substitute equations (3.103) and (3.101) into (3.78) for CP 



case, we obtain 



N-l 



N k = (k + l)(N-k-l) 



(3.104) 
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Chapter 4 

SU(N) Bogomolny Solutions 



4.1 Harmonic Maps Ansatz 

In this chapter, we discuss the general constraction of the solutions for SU(N) Bogo- 
molny equations for BPS magnetic monopoles using previous description in harmonic 
maps. We also give an example for SU(2) Bogomolny equations that give the same 
solutions with one was proved by M.K. Prasad and CM. Sommerfield |17j . 

For this purpose, we take an ansatz for SU(N) Bogomolny equations given by T. 
Ioannidou and P.M. Sutcliffe HHI 



N-2 , y o N ~ 2 



$ = i £ h 3 {r) P, - - , = 9j(r) [PjAPj) > A r = (4-1) 

3=0 V JV/ j=0 

We substitute the ansatz into SU(N) Bogomolny equations (2.47) and (2.48) then 
multiply from the right with vectors field Mi as follows 

(zD e $ - F r ^)Mi = (4.2) 

(*D r * - (1+ Jf )2 %) Mi = (4.3) 

Next, the equation (4.2) is expressed in term of projection operators Pk- In that case, 
we need to describe the properties of projection operators Pk and their derivatives 
applied to Mf. 

P k M t = 5 k M (4.4) 
from equations (3.59) and (3.60), we obtain 

d^PM = 9 € (M fc |M fc |- 2 )A<EM, + M fc |M fc |- 2 ^MjM, 
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= (hi - M i+i 



(4.5) 



d^P k M t = dz(P k M{) - P^Mt 

= (5 k ,i-i - 5 ki ) M^i/Ci-i (4.6) 

where /Q_i = |M/_i| _2 |M/| 2 . For simplicity, we use a convention for index j which 
means the sum of j — 0, . . . , TV — 2. 

Therefore, we can take the action of the fields, as given by the ansatz (4.1), on 
the vectors field M t by using properties (4.4)-(4.6) 

A^Mt = -G t M l+1 (4.7) 

AgM l = G W M W /C W (4.8) 

where Gi — g% + gi-i- Then, we derive that 

F^M t = -d r GM +1 (4.9) 

F a M t = (deAf - <%4 e + [At, Af] ) Mi (4. 10) 

For equation (4.10), we calculate separately for its parts as follows 

dtAfMt = M,(/Cj_iGi_i - Kid) (4.11) 

%4 C M, = M,(/C,G, - K H G H ) (4.12) 

A^M t = -G?_iM,X;,_i (4.13) 

A f ^M, = -GfMiKi (4.14) 

Substitute equations (4.11)-(4.14) into equation (4.10), brings 

= M, (/C,_i Qi-! - /Q Qi) (4. 15) 

where Qi = G t (2 - Gi), then 

§M l = i(h l -^M l (4.16) 

d r QM l = i(^ r h l -^^jM l (4.17) 

9 € $M, = - h l+1 )M l+1 (4.18) 
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Substitute all the previous result into equation (4.2), gives 



{i(d € $ + [^,$])-F rC }M, = 
{(h l+1 - h t ) - (h l+l - hfiGt + d r G t }M l+1 = 

because M^ +1 7^ 0, therefore 

(hi =l -hi)(l-Gi) + d r Gi = 

If we do the same way to the equation (4.3), then we have 

Vi - ^1 M, - (1 tlf ^ M^^-xQ,-! - /C,Q,) = 



(4.19) 



(4.20) 



N 



2r 2 



(4.21) 



For CP( N ^ case, we take vector fields Mi as given by the Veronese map from the 
previous chapter 



i,vcf^,...,vcrr,...,e 



then we may write the factors (1 + |£| 2 ) 2 /Q as constant K h such that 



JC, 



2\2 



(4.22) 



(4.23) 



where we find that Ki = Mi = {I + 1){N — l — l) are constants which depend on index 
I as shown from the equation (3.100). Hence, the equation (4.21) becomes 



' d r hj 



drh^J-^ Q,_i - K x Q t ) I M, = (4.24) 



because M; 7^ 0, then we obtain 



(^p- - d r h^j - ± (K^ Qi^ - Ki Qi) = (4.25) 



or in complete 



Sj=0 drhj 
N 



d r h t 



2r 2 



(Ki^Qi^-KtQ, 



(4.26) 
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Note that, by definition h h g x = if l£ {0, 1, . . . , (N - 2)}. 

The equations (4.20) and (4.26) are called SU (N) Bogomolny equations for BPS 
magnetic monopoles. Clearly, those equations are in of scalar and non-linear coupled 
fields. Compared to the form of the SU (N) Bogomolny equations before we apply 
harmonic maps method, the previous SU(N) Bogomolny equations are in form of 
matrix and aslo non-linear coupled fields. It is much easy to work with scalar field 
equations than in form of matrix field equations. But, to find the exact solutions 
of SU(N) Bogomolny equations, even in scalar field equations, is another complex 
problem since the fields in form of non- linear coupled. In this thesis, we do not discuss 
about the exact solutions since it needs others advance mathematical techniques to 
find the soliton solutions. 



4.2 SU(2) Bogomolny Equation 

In this section, we show a simple example for 577(2) Bogomolny case for which the 
equations (4.20) and (4.26) has range of value I — 0, . . . , (2 — 2) = 0, so by definition 
h h gi — if I 7^ 0, then equations (4.20) and (4.26) become 

- h (l - g ) + go = (4.27) 
-h + ^g (2-g ) = (4.28) 

with g = d r go and h = d r ho 

Now, define a function K — 1 — g and substitute it into equations (4.27) and 
(4.28), such that 

- h K - K = (4.29) 
-h + ^(l~K 2 ) = (4.30) 

From equation (4.30), we have 

h = ±(l-K*) (4.31) 

Take the differential of equation (4.29) about r and substitute it into equation (4.31) 

- h K - h K - K = 
K = — \k(\ — K 2 ) — h K (4.32) 
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from equation (4.29), gives 

r 2 K = KiK 2 - 1) + r 2 h 2 K (4.33) 

Take the differential of equation (4.30) about r and substitute it into equation (4.29), 
then 

- h + \(-2KK) -hl-K)=0 

r 2 h = --(1 -K) + 2h K 2 (4.34) 
r 



Next, we define a function H = tHq with 



h = 


H 




H 




+ 








r 


2H 


2H 


+ 


H 






r 



(4.35) 
(4.36) 



and substitute into equations (4.33) and (4.34) 



r 2 K = K(K 2 - 1) + H 2 K (4.37) 

r 2 H = 2K 2 H (4.38) 

The result in equations (4.37) and (4.38) are same as the one that was obtained by 
M. Prasad and C. Sommerfield p2] which gives 

K = (4 - 39 > 

H = C rcoth(C r) - 1 (4.40) 

where C is a constant. 
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Chapter 5 



Conclusions 



In this thesis, we have described how to use harmonic maps in CP( N ~^ space to 
simplify the SU(N) Bogomolny equations for BPS magnetic monopoles. As the 
result obtained in chapter IV equations (4.20) and (4.26), we just need to solve the 
non-linear coupled equations for scalar fields that depend on variable r, as follows: 



with I = 0,1,..., (TV -2). 

Even it is not a simple problem to solve the non-linear coupled equations, but at 
least we have more simple form of the equations in scalar than before. As an example, 
we also have shown for £77(2) Bogomolny equations and verified the result with the 
one obtained by M. Prasad and C. Sommerfield. 



(/l J+1 -/l,)(l-G?,)+0rG, = O 



(5.1) 




(5.2) 
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Appendix A 



Riemann Sphere 



In deriving metric (2.19), we introduce equation £ = tan(|)e 1 ^ and transform vari- 
ables of the metric (2.18) into complex variables (9, 0) — > (£, £). By using other forms 
of the transformation equation which are 

tan (fU 7^ = K|, sin ( 6 -) = -=JH_, cos f ?■) = , 1 ; = (A.l) 



,so that 



(i+iei 2 )r 5 ^ (i + iei 2 )e 

%>=-^, %/> = ^ (A.2) 
with the result of equations (A.l) and (A.2), then d6 2 and d(j) 2 are written as 



« + ¥ = (1 _^ |2)2 + |j^ 2 J + (1 + |e|2)2 (« + «) 

rf£ + «%0 rf£) 2 = -1 (^£ 2 + Id? - ±^(d£dZ + dtd^ (A.3) 



and 

sin 2 ^ = 4 sin 2 j cos 2 |^ 



4 ^ (A.4) 



(1 + iei 2 ) 2 
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Substitute all the result from (A. 3) and (A. 4) into equations (2.18) then we have 
metric (2.19) 

2r 2 

ds 2 = dt 2 - dr 2 - - rTT ^(df + 



(1 + iei 2 )^ 



So, we can write its tensor metric by 



(A.5) 



9w 



( 1 o 

-1 



V o o 
/ i o 

-1 



V o o 







2r 2 



(i+l?l 2 ) 2 




(1+lgl 2 ) 2 






2r 2 



2r 2 



(i+ICI 2 ) 2 





(l+l£l 2 ) 2 

2r 2 





(A.6) 



(A.7) 
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Appendix B 



Derivation of Equations (3.40) and 
(3.41) 



To derive equation (3.40), we need to use tensor metric (A. 7), such that equation 
(3.39) becomes 

[P,d^P] = [P,g^d^P] (B.l) 

As stated in section 3. .4. 2 that the source manifold is a two dimensional Euclidean 
space which is wrote is Riemann sphere coordinates. It means that tensor metric 
(A. 7) only use variables of Riemann sphere (£,£), then the equation(B.l) becomes 



[P^d^P] = [P,g^d^P + g^d^P] 



(1 + l£l 2 ) 2 



.,2 



[P,d&P] (B.2) 



and from equation (3.39) then 

[P,dftP]=0 (B.3) 
Next, for equation (3.41), we write equation (3.40) as follows: 

[P,dftP] + [P,dtd i P] = 
<%(P<%P) - %(%PP) + <%(P<%P) - <%(<%PP) = 

[P, %P] + %[P, %P] = (B.4) 
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Appendix C 



Derivation of Equations (3.49) and 
(3.50) 



Equation (3.49) is derived by writing it in form of matrix field M: 



and if we take its Hermitian conjugate, then 



While for equation (3.49), we can derive it from instanton solution as follows: 



Pd^P = 

Pdz(I-P) = 

d ( (P(I - P)) - d ( P (I - P) = 

d^{M\M\~ 2 M\l -P))-d^{M\M\- 2 M^){I -P) = 

-d^M\M\~ 2 )M\l -P)-M\M\~ 2 d^M ] {I -P) = 

M\M\~ 2 dt.M\l - P) = 



(I — P)M 



M — M\M\~ 2 M^M 




(C.l) 



((/-P)M) f = 
M\l-P) = 



(C.2) 



M\M\- 2 {d^M) ] {I - P) = 



(C.3) 
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Appendix D 



Derivation of Properties (3.59) and 
(3.60) 

In this Appendix, we present derivation of the properties of equations (3.59) and 
(3.60). For equation (3.59), we use the fact that the sequence of M k are independent, 
so we have the expansion: 

A 

<%M z = ^M fc a fcZ , / = 1,...,A (D.l) 

fc=0 

where A k i are the (n x n) matrices. We have an assumption that M is holomorphic 
<%M = 0. 

Multiplying (D.l) from the left by then we have 

M^-M, = \M m \ 2 a ml (D.2) 
From recurrence relation (3.56), gives 

M m+1 = {I-P m )d x iM m 

Mi +l = (d^M m y(i-p m ) 

(%M m )t = Ml +1 + (^M m )t Pm (D.3) 

and the left side of equation (D.2) becomes 

M^-M, = d^MlMt) - d^MlM t 

= d l {M}M l 8 lm )-{d^M m )^M l 

= 8 lm M}d i M l + 5 lm d- i M}M l -{d i :M m yM l (D.4) 
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Substitute equation (D.3) into (D.4), such that 

M^-M, = 6 lm M}d^M l + S lm d^M}M l -(M} +1 + (d i M^P m yM l 

= 5 lm M}d^M x + 5 lm d^M}M l - 5 l<m+1 \ M t \ 2 - d^M^ 

= S lm Mld^M t + 5 lm d-^M}M l - £, )m+ i|M,| 2 - (J^^-m/m, 

= 5 im M/^-^ m+1 |M z | 2 (D.5) 

and substitute the result into equation (D.2) as 

M^M t = \M m \ 2 a ml 

= 5 lm Mld^M l -5 l , m+l \M l \ 2 
a ml = 5 im |M m |- 2 M/%M, -<5,, m+1 |M m |- 2 |M i | 2 (D.6) 

If we substitute equation (D.6) into (D.l), we obtain 

A 
fc=0 

= £ M fe (<J, fc |M fc |- 2 M;^-M, - ^ +1 |M fc |- 2 |M ; | 2 ) 

fc=0 

= M^M^Mjd-^-M^M^-^Mtf (D.7) 

Because a zz = |M,|- 2 M/%/kf; and aj_i,, = — | A^ z _! | - 2 1 Af z | 2 , then 

= M,a„ + M^a^ (D.8) 

Next, we prove that an = 0, / = 1, . . . , A or equivalent to Mjd^Mi = 0. For the case 
I = 1, using the construction (3.57) or (3.58), we obtain 

M 1 = (I - P )^M , <%M = 

<%M Z = d i (I-P )d^M + (I-P )d^M 

= -<%P <%M 

= -d^M \M \- 2 M^M 

= -M d^(\M \- 2 M^M 

M/<%M z = (D.9) 

For general case with 1 < k < A, we use the inductive argument by assuming that 
a kk = such that <%M fc = M fc _ia fc _i ife . 
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As 

«%P fe = d^M k \M k \~ 2 Ml + M k d-^\M k \- 2 Ml) 
= (d^M k + M k d^)(\M k \- 2 Mt) 

= (M fe _ 1 a fc _ life + M k d^){\M k \- 2 Ml) (D.10) 

using the construction (3.59) and orthogonality property (3.51), it follows that 

M k+l = (J-P P k )dl +1 M 

<%M fc+1 = <%(/-P Pk)dl +1 M + (/ - P P fe )<9f +1 <%M 

= -d^P + --- + P k )dl +1 M 

= -^%P4 +1 M (D.ll) 



1=0 



If we multiply equation (D.ll) from the left by M k+1 

%P4 



1=0 

k 



o 

1=0 

= (D.12) 
then it shows that a k+lik+ i = 0. Therefore, equation (D.8) becomes 

<%M, = M iail + M^cn-u 

= -M^lM^r 2 ^! 2 (D.13) 

which is proof of equation (3.59). 

We may write equation (D.13) as follows: 

d^Ml = K-i^M^ (D.14) 

and use the orthogonality property (3.50), then 

d ( \M k \ 2 = M\d^M k + d^MlM k 

= M^M k + {a k _ hk y Ml^M k 

= M ] k d^M k (D.15) 
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such that 

d^M k \M k \~ 2 ) = d ( M k \M k \~ 2 + M k ds\M k \- 2 

We have 

^(|M fe | 2 |M fc |- 2 ) = 
\M k \%\M k \- 2 + d ( \M k \ 2 \M k \- 2 = 

d ( \M k \- 2 = -\M k \- 2 d ( \M k \ 2 \M k \ 

and substitute it into equation (D.16), so that 

dt(M k \M k \~ 2 ) = d^M k \M k \~ 2 - M k \M k \- 2 d^\M k \ 2 \M k \~ 

= d^M k \M k \- 2 -M k \M k \- 2 Mld^M k \M k \ 

= (I - P k )d^M k \M k \- 2 

= M k+l \M k \~ 2 

which is proof of equation (3.60). 
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Appendix E 

Derivation of Equations (3.82) and 
(3.84) 

Derivation of equation (3.82) works as follows: 

\M k \^d^M k \M k \- 2 = M fc t ^(-M fc _ 1 |M fc _ 1 |- 2 |M fc | 2 )|M fc |- 2 

= (MlM k ^\M k ^\- 2 \M k \ 2 ) \M k \- 2 
+ {d-^M k )^M k . l \M k ^\- 2 \M k \ 2 \M k \- 2 
= -iMfeHM^rX^M^lM^I- 2 
= -\M k \ 2 \M k ^l 2 (E.l) 

While for equation (3.84), we derive as follows: 

Mld^(P k d ( M k )\M k \- 2 = d^(MlP k d ( M k )\M k \- 2 - d^MlP k d ( M k \M k \- 2 

= d^Mld ( M k )\M k \~ 2 - d-^MlM k \M k \- 2 Mld^M k \M k \- 2 

= % (d 6 (MtM k ) - d^M ] k M k ) \M k \- 2 

= d^(MlM k )\M k \- 2 

= dz (d^\M k \ 2 \M k \- 2 ) - d^MlM k )d^\M k \- 2 

= % (d^\M k \ 2 \M k \- 2 ) - d^MlM k d^\M k \~ 2 - Mld^M k d^\M k \- 2 

= ^(%|M fc | 2 |M fc |- 2 ) (E.2) 
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